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Abstract 

0^ I A sample of hundreds of simulated galaxy clusters is used to study the statistical 

' properties of galaxy cluster formation. Individual assembly histories are discussed, 

■ the degree of virialization is demonstrated and various commonly used formation 
^ I times are measured and inter-compared. In addition, the fraction of clusters which 

have "recently" undergone a major merger or significant mass jump is calculated 

■ as a function of lookback time and interval. The fraction of three- and four-body 

^5 ! mergers is also studied. 
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2 ' 1 Introduction 

Due to their immense size, galaxy clusters can be easily identified with their 
^ ' dark matter halos, and their clustering and number counts can be reliably pre- 

■ dieted by dark matter simulations. However in order to find or "weigh" galaxy 
clusters observationally, assumptions of dynamical equilibrium and less under- 
stood astrophysics usually come into play ^ . Since halos form via the accretion 
and mergers of smaller units, and for clusters this process is occurring to the 
present time, such assumptions need to be specified precisely. For example, 
how well a galaxy cluster is described by equilibrium properties depends upon 
whether it has recently undergone a major merger (where "major" means dis- 
rupting equilibrium) and upon the specific methods used for its detection and 
mass measurement. Thus even if one is only interested in galaxy clusters as 
"test particles" or peaks in the density field, the history of cluster growth is 
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important in order to make contact with observations. Cluster growth histo- 
ries are also important for understanding other cluster properties. For exam- 
ple, several observational phenomena (see below) are associated with clusters 
which have recently undergone mergers. A related question is which cluster 
properties have the least sensitivity to cluster assembly histories. 

In this paper we study the assembly history and degree of virialization of high- 
mass halos in large N-body simulations. A cluster assembly history can be 
characterized either by events with specific occurrence times, such as mergers 
or large mass changes, or by properties of its entire history, e.g. a parameteriza- 
tion of the mass as a function of time. We calculate several such quantities for 
a statistically significant sample of halos. We consider several popular forma- 
tion time definitions and cluster history parameterizations and statistics. We 
also calculate the fraction of "recently" merged galaxy clusters as a function 
of redshift, back to z ~ 1. 

The definitions of both "merger" and "recently" depend on the cluster prop- 
erty of interest and we explore several choices. We also find the fraction of 
galaxy clusters which have had a recent large mass gain (including accretion) 
for several choices of interval and two final/initial mass ratios. Merger histories 
can be reliably extracted from N-body simulations, and as such these recently 
merged fractions are implicit in earlier work. However, it can be difficult to 
obtain specific numbers from the literature for ACDM models, especially if 
one has a particular relaxation time in mind. In part this is because previous 
studies of different quantities at different times have been published over sev- 
eral years, often using different cosmological models. Here we compute several 
of these quantities for a much larger sample than used in earlier work, for 
ACDM cosmologies, and present them in a homogeneous manner in the hope 
that this will be a useful reference for the community. 

The outline of the paper is as follows. Section 2 is a review, including pointers 
to earlier work on formation times and examples of observed merger phenom- 
ena. Section 3 describes the simulations and methods. A reader interested 
primary in the results can skip directly to Section 4, which has comparisons 
and distributions of some cluster formation properties, and the fractions of 
clusters which have recently merged or had a large mass increase, as a func- 
tion of time. Three and four-body major mergers are also studied. Finally, 
Section 5 presents our conclusions. 



2 Background 

The growth of structure by mergers and accretion is key to the hierarchi- 
cal paradigm of structure formation, and thus mergers and mass gains have 
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been studied intensively. Previous work on cluster formation histories includes 
Refs. (2; 3; 4; 5; 6; 7; 8; 9). While cluster assembly is a complex and ongoing 
process in hierarchical models, it is often useful to have some measure of when 
the cluster "formed". Refs. (5; 7; 8), each with 10-20 clusters, considered for- 
mation time definitions including the redshift, -Zjump, of the most recent large 
AM over a short time. Both (7; 8) also found a characteristic formation time 
Zf associated with the entire cluster growth curve using the parameteriza- 
tion of Ref. (10). This parameterization works extremely well for galaxy sized 
halos and correlates with other properties such as concentration. For galaxy 
clusters, Ref. (8) introduced a generahzation to help better match the more 
recently active formation histories of galaxy clusters. Ref. (6) had a sample 
similar to ours and found a "turning point" time where halos went from a 
quickly growing phase to a more slowly growing phase, this turning point was 
correlated with concentration. There are other formation times considered in 
the hterature. For instance Ref. (11) found the average mass accretion history 
of halos generated by extended Press-Schechter. Ref. (7) also looked at the 
amount of mass gain coming from large AM "jumps" as fraction of cluster 
mass. We consider these properties and their distribution for our large sample 
of over 500 clusters for a ACDM erg = 0.8 model in §4. 

Our second set of results is the frequency of recent mergers and recent large 
mass gains for some fixed lookback time and definition of "recent." Analytic 
estimates of merger rates, in many cases combined with simulations, include 
those by Refs. (12; 3; 5; 13) (see also the extension reported in (14)). The 
previous work closest to the merger counts considered here is by Ref. (5) who 
found the major merger distribution as a function of redshift for a system 
which had 11 clusters aX z = (with a merger defined as an increase in halo 
mass such that Mf/M^ > 1.33) and that of Ref. (13) which also studied the 
recently merged population, but for a smaller number of clusters and not in as 
much detail. It should be noted that merger rates and the number of clusters 
which have recently undergone a merger are slightly different, it is the latter 
we study here. 

Mergers can alter X-ray temperatures, luminosities, cluster galaxy light, clus- 
ter galaxy velocities and Sunyaev Zel'dovich (SZ; Ref. (15)) "flux" ^ . Ob- 
servational consequences of merging clusters of galaxies have been studied 
analytically and with numerical simulations, e.g. Refs. (20; 21; 7). In these, 
the time scale for major disturbances of clusters due to mergers seems only 
weakly dependent upon the impact parameter (the magnitude of the distur- 



^ Several simulations with different heating and feedback prescriptions have found 
small merger induced scatter in the integrated SZ flux relation (Refs. (16; 17; 18) 
although there are examples where a ^ 50% increase in total SZ flux for one sound 
crossing has been seen (19)). Significant scatter in the mass-flux relation instead 
appears to be dominated by that due to line of sight projection (16). 
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bance has a stronger dependence, e.g. sec Refs (21)). We do not consider 
impact parameters below. The errors induced by ignoring merger effects (on 
X-ray luminosities and temperatures) for cosmological parameter estimates 
has been considered in (22) ^ . More generally, unaccounted for mergers can 
disguise cluster masses and alter survey selection functions. 

Observationally, some signatures of mergers can be used to flag unrelaxed 
clusters. For instance substructure, the presence of multiple peaks in the clus- 
ter surface density on scales larger than the constituent galaxies, has been 
studied extensively. Several different observational methods have been used 
to quantify the amount of substructure in X-ray clusters (see the review (24) 
and references therein). The fraction of clusters exhibiting substructure has 
been found to be significant in several surveys. For example, Ref. (25) visually 
found substructure in 41% of 208 Einstein IPC images, Ref. (26) measured 
the emission weighted centroid variation for 65 Einstein clusters to get a sub- 
structure fraction of 61% (see also (27) for related work) and Ref. (28) found, 
using three different indicators, substructure in 52 ± 7% of 470 clusters in 
the ROSAT all-sky survey. Substructure has also been searched for in other 
wave-bands. For example, substructure was detected in over 80% of 25 low 
richness 2dFGRS clusters ((29); see their references for other studies). Fol- 
lowing Ref. (30), substructure has been considered as a means of constraining 
cosmology, e.g. in (31) using 40 Chandra clusters over a range of redshifts (see 
also (32) for relevant simulations). A meaningful comparison of our cak:ulatcd 
number of "recently merged'' clusters to the number of observed galaxy clus- 
ters with substructure depends on relaxation times. These times depend upon 
the substructure measurement in question and are still not well estimated. For 
instance, relatively small sub-clumps may take a different amount of time to 
disappear depending on their density, see (24) for further discussion. 

Radio halos and radio relics have been proposed as another merger indicator 
(see the reviews (33) and (34)). Clusters hosting radio halos and rehcs tend 
to show other signs of merger activity, however not all clusters with other 
indicators of recent mergers show these radio signals. In a complete sample (35) 
at the NRAO VLA Sky Survey surface brightness limit, only 5% of clusters had 
a radio halo source and 6% had a peripheral relic source, however in clusters 
with X-ray luminosity larger than 10^^ erg s~^, 35% had radio halos or rehcs. 
In the ROSAT sample considered above by Ref. (28), 53 clusters of the 470 
had radio halos or relics, most of these were high luminosity {L > 4.0 x 10*^^ 
erg s~^). It has been proposed that more radio halos and relics will be found 
in lower mass clusters as sensitivities improve (e.g. Ref. (19)). The relaxation 
time for radio halos and relics is not known. 



^ They used an analytic model (23) to get a merger history, in principle our calcu- 
lations here could also be used for this purpose. 
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The absence of mergers (i.e. a relaxed cluster) has been suggested as a re- 
quirement for cool core clusters (e.g. Refs. (36)). The fraction of cool core 
clusters has been measured for various samples, e.g. Ref.(37) found with Chan- 
dra that 55% of 38 clusters in the ROSAT Brightest Cluster Sample (with 
2; ~ 0.15 — 0.4) show mild coohng (associated with a coohng time less that 
10 Gyr) and 34 % show signs of strong cooling (a cooling time less than 2 Gyr). 
There are some simulations that suggest however that mergers may help form 
cool galaxy cluster cores (e.g. see (38)), also 5/22 of the ROSAT cooling flow 
clusters have significant substructure (28). 

There are many more phenomena associated with mergers which are even 
more complex (see Sarazin's list of questions about merging clusters of galax- 
ies (39)). The examples of merger indicators reported above have had their 
occurrence rates measured in some of the largest cluster samples available. 

There are some differences between mergers and any form of mass gain ^ . 
Clusters which have gained mass primarily via accretion might be expected 
to be more relaxed than those who reach the same mass via a major merger. 
For smaller halos those with a relatively large mass gain over a short time are 
more clustered than halos of the same mass but without this mass gain (41). 
The subset of recently merged halos, in contrast, does not appear to be so 
biased (42). The authors of Ref. (41) conjecture that objects in denser regions 
tend to have more nearby material to accumulate, resulting in a bias, while 
the same is not true of recent mergers. 



3 Methods 



To investigate these questions we used two N-body simulations run with a 
T^eePM code (43). Each simulation evolved 512^ particles in a periodic cubi- 
cal box, 300/i~^Mpc on a side, for a particle mass of 1.67 x IO^^/i^^Mq. The 
models had h — 0.7, Qm = 0.3, fl\ = 0.7 and erg = 0.8, 1.0. The two normal- 
izations span the observationally preferred range. Outputs were dumped at 
equal intervals^ of conformal time, St = 100/i~^Mpc (comoving), starting at 
z — 2. Other parameters and details of the simulations can be found in (44). 

For each output groups are identified with a friends of friends algorithm (FoF; 
(45)) using a linking length 6 = 0.15 in units of the mean inter-particle spacing. 

^ A modified version of the PS formalism that differentiates between (instanta- 
neous) major mergers and accretion has been developed in Rcfs. (40). 

This is the same time interval considered by Rcf. (7) and shorter than that con- 
sidered by Refs. (5; 8), though the latter groups stiU had their time spacings weh 
within the lower limits for a merging time scale. 
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These groups correspond roughly to all particles above a density of 3/(27r6^) ~ 
140 times the background density. We did not consider other group finders, a 
comparison between HOP (46) and FoF merger fractions for a smaller sample 
was considered in (13). In that case there were both fewer recently merged 
clusters and fewer clusters altogether for the FoF case. See also Ref. (41) for 
a comparison of HOP and FoF mergers and accretion for smaller mass halos. 
They estimate what minimum mass fraction in accretion is robust for different 
group fractions and find that a 20% mass increase (for a shorter length of time 
and smaller halos than we consider) is outside the "noise" region for the FoF 
vs. HOP group finders. 

For each group and output time we compute a number of properties, including 
the mass, velocity dispersion and potential energy of the group. All properties 
are computed using only group members. Masses henceforth refer to friends 
of friends masses Mfof unless otherwise stated. To get M200 or Myi^ conversion 
is needed (47) - for a fitting formula see (48). From the simulation we find 
M = 1.2 M200 at 2 = 0, and M = 0.98 M200 at 2 ~ 1. We consider halos to be 
clusters if M > 10^^ /^"^M© (i.e. 5971 or more particles); at ^ = there are 
909 clusters in the as — 1.0 box and 574 in the as = 0.8 box. 

For the formation time studies, most of the definitions are straightforwardly 
derived from the cluster histories. We track the properties of all of the pro- 
genitor halos with more than 64 particles for all of our clusters. From these 
histories it is straightforward to find the mass of the largest progenitor as a 
function of time. To find mergers, the two largest predecessors are found for 
each cluster at an earlier time. The final cluster is said to have had a (major) 
merger between the earlier time and the later time if the ratio of particle num- 
bers going from both of the predecessors into the final cluster is above some 
minimum ^ . In the case where the time interval of interest does not coincide 
with a dump time for the data, the two closest times bracketing the time of 
interest are taken, and the merger fraction is linearly interpolated between 
these (starting from the time closest to the time of interest). One relaxation 
time we consider is the crossing time, in this case we start at the earlier time 
and interpolate between final times. 

We checked our definitions of mergers in two ways. Instead of using the mass 
of only progenitor particles which actually go into the final halo in the mass 
ratio, one can use the full mass of the progenitor halos, this changes the merged 
fraction by less than 5% for the longer time intervals and almost always less 
than 10% for the short time intervals. We prefer our definition as it might 
help counteract the inclusion of many unbound particles in halos due to our 
FoF group finder. Secondly, we identify mergers by looking at the initial and 



This was in part to try to minimize the effect of particles which were not actually 
bound being taken to be part of a cluster's predecessor. 
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final halos and comparing mass ratios of the initial halos, without tracking 
their intermediate behavior until they merged to the final halo. A concern was 
that two predecessors with an initial mass ratio of e.g. 1:5 might first become 
a mass ratio of 1:10, say, if the largest progenitor grew significantly in mass 
relative to the second largest, before merging. Looking at intermediate times, 
the mass ratio of predecessors is differs for the two definitions, but mostly for 
almost equal mass predecessors. Once the predecessor mass ratio drops to ~ 
1:3 the scatter between the two ways of defining predecessor mass ratio gets 
quite small. Using the more accurate definition based on intermediate time 
steps produces a slightly larger number of recently merged clusters, again less 
than 5%. That is, our approximation in some cases underestimates the number 
of mergers by a small amount. 

We vary the mass ratios and time intervals, and consider intervals equally 
spaced in real time, in light crossings (scaling as a) and in terms of "cross- 
ing times" (scaling as (Gp)^^/^ ~ o;^/^). We call the lookback time the time 
at which the cluster is observed, and the relaxation time the interval prior 
to the lookback time within which the merger or mass jump ("large AM") 
has or hasn't occurred. Intervals of 100, 300 and 600a/i~^Mpc hght cross- 
ings correspond to approximately 1, 3 and 6 x a/2Gyr. We call these St — 
100, 300, 600 /i~^Mpc in the following. Although the simulations go back to 
z ~ 2, we only plot back to lookback times where the statistics have any 
constraining power. 

We also consider three- and four-body merger fractions, since such systems 
have been observed (see, e.g. Refs. (49)). We define three-body mergers as the 
case when the second and third largest predecessors each contribute to the final 
halo at least 20% of the mass contributed by the largest predecessor. Similarly, 
four-body mergers require the fourth largest predecessor to also contribute at 
least 20% of the mass contributed by the largest predecessor. 



4 Results 

4-1 Formation times and properties 

A classic paradigm for cluster assembly histories is the spherical top-hat col- 
lapse model. In this model a uniform, overdense region "breaks away" from 
the expansion of the universe, evolving as a self-contained positive curvature 
universe. This toy model is often used to motivate the threshold density for 
virialization or to estimate formation times (50). The evolution of such per- 
turbations does not closely resemble the rich structure seen in the formation 
of massive halos in N-body simulations, as noted by many authors referenced 
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above, however it is interesting to see how it performs nonetheless. We show 
in Fig. 1 the evolution of the peak circular velocity (M oc v^) of the most 
massive progenitor for 6 clusters, chosen at random from amongst the most 
massive (M > 6 x 10^^/i~^Mq) in our — 0.8 simulation. Since the mass in 
the spherical top-hat model is difficult to compare to simulations (it is con- 
stant and has uniform overdensity) , we use peak circular velocity. The peak of 
the circular velocity curve, v1 = GM{< r)/r, is computed for each progenitor 
using the minimum of the cluster potential as a center for defining r. The 
dashed line shows the evolution predicted using the spherical top-hat collapse 
model for a cluster which viriahzes at the present. We see that the spherical 
infall description has a cluster's peak circular velocity growing most quickly 
right before virialization, while clusters forming in a cosmological simulation 
have more of a steady growth over time ^ . We have also marked on the plot 
four of the formation times defined and discussed below, 2;juinp, Zi/2-i Zj and 

A different view of the formation of clusters is given in Fig. 2 which focuses on 
the same 6 clusters as Fig. 1. In the top panels we show the mass of the largest 
progenitor as a function of time. In the middle panel we show the ratio of the 
ID velocity dispersion to the peak circular velocity - this has been used as 
a proxy for the degree of virialization of the clusters (51). In the lower panel 
we show the relation between the total kinetic and potential energy explicitly. 
Both energies have been computed using only the particles in the FoF group, 
but no proxy has been used, unlike the middle panel. 

Even focusing only on the most massive progenitor, rather than the full merger 
tree, the upper panels show the rich phenomenology of cluster formation. 
There are several obvious mass jumps interspersed between periods of rela- 
tively smooth accretion. The lower panels show that the cluster is disturbed 
by mergers and large mass jumps with strong departures from the expected 
virial relation 2KE ~ PE. The departures last for several of our output times. 
The "standard" vacuum virial relation 2KE = PE is not satisfied by these 
clusters even in their quiescent phase. The ratio 2KE/PE is typically larger 
than unity due to the steady accretion of material onto the cluster ^° (52; 51) 
and there is an overall decline in 2KE/PE with time. Averaging over the entire 
cluster sample we find for objects of fixed mass the KE/PE ratio decreases 
with decreasing redshift. For objects at fixed redshift the ratio increases with 
increasing mass. For objects at fixed number density the ratio decreases with 
decreasing redshift. For fixed mass cut or number density the decline is steeper 



^ For the spherical top-hat collapse model, the cluster mass is constant and the den- 
sity/radius change with time. In physical coordinates the perturbation first expands 
(falling and then collapses (rising v^. 

This can be thought of as a surface pressure term which modifies the virial rela- 
tion, boosting the kinetic energy. 
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Fig. 1. The evolution of the peak circular velocity of the most massive progenitor 
for 6 massive clusters (symbols connected by solid lines) compared to spherical 
top-hat collapse of a cluster virializing at z = (dashed). The cosmological model 
is our ACDM as = 0.8 model. The large symbols indicate various formation times 
discussed in the text, as indicated by the legend. Only the cluster denoted by the 
filled circle has a positive Zf. 

below z = 1 than above z = 1 and within errors the decline below z = 1 has 
the same shape for all samples. We hypothesize that this decline in KE/PE 
as due to the cessation of structure formation due to A domination below 
z ~ 1 (plus a small shift in the mass of objects currently undergoing rapid 
accretion). We will see a similar drop in merger activity below z ~ 1 in the 
next section. 

We attempted to further quantify the relaxation of clusters to the background 
virial relation as a function of time since last major disturbance (c.f. Ref. (7)), 
but we found the scatter was too large for us to draw robust conclusions. 
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Fig. 2. A view of cluster formation for the same 6 massive clusters as in Fig. 1 (see 
inset legend for correspondence). The left and right columns show 3 clusters each. 
The top panel shows the mass accretion history since -z = 2, the middle panel shows 
the ratio of velocity dispersion to circular velocity (which has been used previously 
as a proxy for degree of virialization) while the lower panels show the virial ratio, 
2KE/PE, as a function of time (sec text). In the lower right panel, the heavy solid 
line is the virial ratio averaged over all clusters with M > 10^^ hT^MQ at the given 
redshift. The first point, z ~ 1.6, has 3 clusters above the threshold. 



Even though cluster formation is a complex process it is sometimes useful to 
attempt to describe it with a single "formation time" . Any such compression 
of information must be imperfect, and to some extent arbitrary. Depending on 
the phenomenon of interest different times may be more or less appropriate. 
For this reason several definitions in the literature exist for the "formation 
time" of a cluster. We calculate several of these and their distributions for the 
(78 = 0.8 sample of 574 clusters below. 
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The first definition we consider is when a cluster had its most recent large 
AM, i.e. Mf /Mi > 1.2 in an interval 6t = 100 h~^Mpc. In this section a mass 
"jump" refers to a mass gain of at least this much. For any cluster Zj-amp is the 
redshift of the most recent such jump. In our sample 13 clusters had no large 
mass jumps at any time after z ~ 2, two of these clusters are more massive 
than 3 x IO^^/i-^Mq. Ref. (7) found that 2 out of their sample of 20 did not 
have any large mass jumps. Another common definition is when the cluster 
has reached at least half its mass, for this time there is an analytic formula 
(53) which was tested in (and improved using) simulations in Ref. (54). We 
define Zi/2 as the earliest output time when the cluster had at least half of 
its mass. The relation between these two definitions is shown in the upper 
left hand panel of Fig. 3. There is a large amount of scatter between the two 
definitions but there is a striking transition at z = 0.5. For clusters which 
have not reached half of their present day mass before z — 0.5, almost all have 
had a large AM jump after reaching half-mass, i.e. they have at least half of 
their mass before they have their last large mass jump. However, for clusters 
which attained Zi/2 before redshift z = 0.5 (286 clusters), slightly under one 
third (80) had their most recent mass jump even earlier (i.e. ^ijump > ^1/2)- 
An intuitive reason for this is that clusters which reached half of their mass 
early on have not been gaining mass quickly and so are less likely to have 
large mass jumps after zi/2. This general shape persists if one considers some 
other "formation time" such as 2:1/3 or Z3/4 (with obvious definitions), the z 
value where -Zjump starts becoming larger (that is, earlier) than 2;i/3(2;3/4) is 
larger (smaller) than that for zi/2- A larger fraction (12/33) of the clusters 
with mass greater than 3 x 10^^ h~^MQ and Zi/2 > 0.5 had ^jump > ^1/2- H 
one is interested in formation times in order to find undisturbed clusters, for 
lower redshifts it appears that -Zjump may be a more conservative estimate. 

A third definition of a cluster formation time is when the potential well of 
the object becomes deep enough to be considered "a cluster" . We take this, 
somewhat arbitrarily, to be when the object has reached 10^^ /i^^Mq, and call 
this zu. 

All of these definitions rely on special events in the mass accretion history. An 
alternative is to use the whole history of the most massive progenitor. While 
this is less information than in the entire tree, it provides a global view of the 
formation process. We consider three fits used in the past for simulations. For 
galaxy sized halos, Ref. (10) found that the mass accretion histories could be fit 
by M(a) = Mqc"^"^^ = Moe~^^/(^+^^\ and that Zf correlated well with other 
cluster properties, i.e. concentration. Ref. (8) found that for clusters they could 
obtain better fits if they generalized the function to M(a) = MoaPe~^^f^ ~ 
dPM{a)^ where aj = 1/(1 + 5/). A third form was proposed by (6), who fit the 
the mass accretion history to the form Mj/Mtp = / (pvir(-2tp)/Pvir(-2t)) where 
Pvir^z) is the virial density at redshift z computed from spherical top-hat 
collapse and / interpolates between two power- laws (6). The latter authors 
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found that halos generally transition from a period of rapid accretion to a 
phase of slow accretion and used the transition redshift, Ztp as a proxy for 
"formation" time. 




Fig. 3. Top left: -^^jump- The 13 clusters with no large AM in the simu- 

lations are not shown as they have Zjump > 2. Top right: zi/2 vs. Zf. Lower left: 
Ztp vs. Zjump- Lower right: ztp vs. Zf. The red circles correspond to clusters with 
M > 3 X 10^^ h~^MQ. The correlation coefficient, r, for each case is shown at upper 
left. 

We fit these three parameterizations to our sample, estimating Zf, Zf, Ztp through 
a least squares fit of ln(Mj/Mo) against Zj, with all points being equally 
weighted. The values of Zf, Zf are allowed free range, however Ztp, due to 
its definition, is constrained to the range < z < 2 where we have data . 
We compare this formation time Zf with 2:1/2 at upper right in Fig. 3; they are 

Also, we fit M not M^[j. but this should have little effect since M^i,. ~ 1.08 Af at 
z = and ~ 1.10 M at high z. 
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correlated (as also found by Rcf. (8)). We show in Fig. 3, bottom, the com- 
parison of this "turning point" time z^p with ^jump and Zf above. The largest 
correlation shown is between Ztp and Zf, as shown in Fig. 3. 

Although these parameterizations do describe general trends in the cluster 
mass histories, good fits were not found for all the cluster histories. (The best 
fits using our least squares criterion was for the models of Ref. (8), which 
have an extra parameter. However, similarly to (8), we found that many of 
the clusters had either p = (278/574) or a/ = (119/574).) This made 
results sensitive to the z-range used and the interpretation of the fit coefficient 
difficult. To illustrate the difficultly, we show in Fig. 4 examples of fits for 3 
clusters. 

In more detail, for the fit for Zf (introduced by Ref. (10)), we investigated the 
dependence on the range of z used in the fit. Ref. (8) fit their Zf to histories 

for < z < 0.5 and for < 2; < 10 and found that the results were consistent. 
We found extremely large variations comparing < z < 0.5 and < z < 2.07. 
Our fits may have differed in success from theirs for two reasons. Our clusters 
tend to be more massive (theirs were between 5.8 x lQ^^h~^MQ < Migopj, < 
2.5 X IO^^/i^^Mq, which are even lower masses in terms of our mass definition), 
as they note, less massive halos tend to have better fits in general. Secondly, 
our sample is a lot larger, allowing us to sample a broader range of behavior. 
Even some clusters which gained more than 90% of their mass in the range we 
consider have Zf which depends sensitively on the fit range. This is a simple 
consequence of the fact, which can be seen above, that the functional form does 
not well describe the shapes of the mass accretion histories of these objects. 
Finally we note that for clusters with recent mass jumps the fit prefers Zf < 0. 
As it is very plausible to say that these systems arc still in the process of 
forming we report these values above without renormalization. 

The other two parameterizations also had difficulties with some subset of the 
clusters. It should be noted however that our redshift range was smaller than 
(8) and that they did not report robustness of the p 7^ fits to smaller z 
ranges. Ref. (6) actually considers a sample of similar size to ours for their 
fit, and large scatter is evident around the parameterized fit in their results as 
well, they find a useful correlation between Ztp and concentration nonetheless. 

In addition to comparing the different formation times to each other, we also 
use the cluster sample to get a snapshot of the formation time distribution. 
In Fig. 5 we show the distributions of these different formation times for all 
the clusters at 2; = 0. Clusters which had no large mass jump since z — 2 are 
put in the -Zjump = 2 bin for completeness. 

There is quite a spread in these formation times, however some trends are 
clear. By 2; = 0.2 over 85% of the clusters have had their last large mass jump 
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Fig. 4. Examples of parameterized fits to three clusters, one of the better fits (top), 
less good (middle) and worst (bottom). The solid black line in each panel is the 
history of the most massive progenitor, the solid blue/dotted/red lines are the fits 
of Refs. (10), (8) and (6) respectively. For the top cluster the fits of Refs. (10) and 
(8) coincide. The stars are the associated formation times (for the bottom cluster 
the formation time of (10) is in the future). We list the percentiles in goodness of 
fit for each of the fits at upper left. 

and over 81% of them have formed according to the definition of Ref. (10), 
while 88% of them have at least half of their mass. The different definitions 
have different biases though, e.g., of the redshift formation times shown, only 
the definition of Ref. (10) can go negative. 

A major part of our difficulty in fitting smooth curves to the mass accretion 
histories is the presence of significant mass jumps. These jumps are evident in 
the examples and are a general feature in cluster formation. The statistics of 
these jumps are of interest. Fig. 6 gives the distribution of the number of times 
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Fig. 5. Distributions of different formation time definitions for all clusters in the 
(Tg = 0.8 sample. On left are the distributions for zjump (solid line), 2:1/2 (dotted line) 
and (dashed line). At right, is the distribution for zj (solid line), ztp (dotted 
line) and zj (dashed line). Note that some clusters have zj < 0. See text for further 
definitions of these formation times. 

since z = 2 that the mass of each cluster changes by 20%. This is the same 
quantity considered by Ref. (7) for their sample of 20 clusters more massive 
than 1.2 x 10^^ /i~^Mq - our sample is more than an order of magnitude 
larger allowing us to more completely characterize the distribution. We find 
that there is little dependence on the cluster mass - lines for clusters above 
1, 2 and 3 x 10^^ h~^MQ are shown. The average (median) number of jumps 
for each of these samples is 4.0 (4), 4.3 (4), 4.6 (5), respectively, with a wide 
spread. The right panel of Fig. 6 shows how much of the cluster mass is gained 
in these events - slightly less than half of the clusters (272/574) get at least 
half of their mass in jumps of 20% of larger. For the M > 3 x 10^^h~^MQ 
subsample the fraction is larger: 41/79. 

With a several large number of jumps common, it appears difficult to find 
good smooth and simple parameterized fits for all of the histories. However, 
given the correlations between the different fit "formation times" and between 
these and e. g. concentration, these fits may be very useful approximations in 
a large number of cases. 



4-2 Time dependence of recently merged or large AM cluster fraction 



The above shows that the present day cluster population has a very wide 
spread of formation times and frequent major disruptions. For observational 
purposes, another question is: how many clusters at a given redshift ( "lookback 



15 



Q r4.r.-;T". , ■ I , ■ ■ I ;- lTriMr±= 

2 4 6 8 10 
number of 20% jumps since z ~ 2 



fraction of mass in jumps > 20% 







0.2 



0.4 



0.6 



0.8 



Fig. 6. Left: Number of clusters with N mass jumps Mj/Mi > 1.2 since z = 2. The 
three Hues are for clusters with with M > 1, 2, and 3 x lO^^/i^^M©. The aver- 
age number (median number) of jumps are 4. 0(4), 4. 3(4) ,4. 6(5) respectively. Right: 
number of clusters with a given fraction of the total cluster mass obtained in mass 
jumps Mf/Mi > 1.2. 

time" ) have had a large disruption recently (within a given "relaxation time" ) . 

The fraction of clusters with a large mass change AM in a given interval is 
shown in Fig. 7, top, as the dotted lines. This is the total mass change, via 
mergers or accretion, within the 1.0 and 2.5 Gyr time intervals. The two lines 
show the fraction of clusters with final/initial mass > 1.2 and 1.33 (these 
ratios were chosen by Refs. (7) and (5)), with the lower line for the higher 
mass ratio. 

We may be interested in the most disruptive mass gains, which occur when 
the progenitor ratios are smallest. We refer to these events as major mergers. 
Figure 7, top, also shows the fraction of recently merged clusters as a function 
of time for minimum mass ratios 1:3, 1:5 and 1:10. For a lookback time of 7 Gyr 
(z ~ 0.8) for instance and the 2.5 Gyr relaxation time, this means considering 
predecessors 9.5 Gyr ago (at z ~ 1.6). 

The error bars are calculated using binomial statistics. If M of the N clusters 
have merged the most likely fraction is f — M/N with variance: 



We use symmetric error bars. One can see here that the number of clusters 



The mean value is (M + l)/(iV + 2). 




M{N -M) + l + N 
(iV + 2)2(Ar + 3) 



(1) 
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lookback time (Gyrs) lookback time (Gyrs) 



Fig. 7. Top: Merger fractions in the previous 1.0 Gyr, 2.5 Gyr, 6t = 300/i~ Mpc and 
^cross) as a function of lookback time to clusters, for ag = 0.8. The solid lines give the 
fraction of clusters (M > 10^^ /i^^M©) whose two largest predecessors contribute 
numbers of particles with ratios of 10:1 or smaller, 5:1 or smaller and 3:1 or smaller 
(top to bottom) at given lookback times. The dotted lines are the fraction of clusters 
which have Mj/Mj > 1.2 and 1.33. For clarity (binomial) error bars are shown for 
only some of the lines. The others are similar. 

which have had a recent major merger increases in the past, reaching a dra- 
matic 80% for 1:10 or smaller mergers 7 Gyr ago (around z ~ 0.83) for a 
relaxation time of 2.5 Gyr. For shorter relaxation times fewer mergers have 
occurred, as expected. For fixed relaxation times the amount of accretion rel- 
ative to merging changes with lookback time. 

In the lower two panels different relaxation times are considered. For some 
phenomena associated with mergers the relaxation times of interest depend 
on the time of observation. We saw in the previous section that departures 
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from the virial relation lasted several hundred Mpc of conformal time. This is 
not unexpected. A typical relaxation time is likely some fraction of the halo 
dynamical time. Since clusters have p ~ lO^Pcrit and characteristic times scale 
as typical timescales should be 0.1 in where tn = is the Hubble 

time. Ref. (7) noted that X-ray disturbances in their gas simulations lasted 
5t — 300/i~^Mpc (scaling as a rather than the crossing time ^ a;^/^). This 
corresponds to roughly 1.3/i~^Gyr at the present, and shorter times in the 
past. Simulations reported in (21) have major merger related disturbances 
lasting approximately one crossing time. Prom (55), we take this to be 1 Gyr 
at the present. These two relaxation times are used in Fig. 7 bottom, with 
that found by Ref. (7) at left and a^/^ Gyr at right. For the latter case, the a 
used is for the initial time, before the merger, as clusters which merge later 
have a larger a and thus longer relaxation time, and thus will also still be 
unrelaxed at the given lookback time. The ratio between the fractions with 
mergers and large AM does not change as much with lookback time compared 
to the case where the relaxation time was independent of the lookback time. 
This suggests that equating large AM jumps with mergers is more reliable for 
phenomena whose relaxation times scale with a. 

In Fig. 8 we compare different samples of clusters, changing cosmology (erg 
1.0) and mass (taking M > 3.0 x 10^'' H'^Mq). At the top are the fractions of 
recently merged or recent large AM clusters for both cig = 0.8 and ag = 1.0, 
for the relaxation time is 6t = 600 h~^Mpc. At present this lookback time is 
slightly over 2.5 Gyr. The error bars are smaller for ag = 1.0 because there 
are nearly twice as many clusters in the sample. The (7g = 1.0 clusters have 
had fewer mergers and mass jumps than those for erg = 0.8; clustering is less 
evolved for dg = 0.8 so a cluster at fixed mass is more likely to be forming 
in our rcdshift range. At the bottom, the recently merged or recent large 
AM clusters are shown for more massive clusters, M > 3.0 x 10^^ h~^MQ, 
with (7g = 0.8, for relaxation times of 5t — 300/i~^Mpc (left) and for 2.5 Gyr 
(right). As there are significantly fewer clusters (79 rather than 574) at z = 0, 
the error bars are a lot larger and are only shown for the top line. Comparing 
with the same plots in Fig. 7, it is seen that more of the massive clusters have 
had recent major mergers. 

Values for ag = 1 are reported in Tables 1-4 for four lookback times, along 
with A^cius so that the errors can be estimated. 

For all relaxation times there is an increase in the fraction of recently merged 
or recent large AM as the lookback time increases, i.e. there were more ma- 
jor mergers and large mass jumps in the past. The fraction depends on the 
relaxation time chosen. An overview is provided in Fig. 9 where we compare 
the ratios of merged and large AM samples near the present to that of about 
6 Gyr ago, for several different relaxation times in our two simulations. A dou- 
bling of the recently merged or recent large AM fraction between the present 
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Table 1 

The number of mergers or large AM events within the previous 1 Gyr as a function 
of lookback time for our as = 1.0 simulation, comparable to Fig. 7 upper left for 
as = 0.8. 



lookback time 


1:3 


1:5 


1:10 


20% 


25% 


^clus 


0.0 


0.18 


0.24 


0.34 


0.59 


0.34 


909 


2.2 


0.22 


0.31 


0.45 


0.73 


0.49 


739 


4.0 


0.29 


0.38 


0.57 


0.83 


0.67 


573 


6.2 


0.45 


0.54 


0.72 


0.97 


0.90 


346 



Table 2 

The number of mergers or large AM events within the previous 2.5 Gyr as a function 
of lookback time for our cts = 1.0 simulation, comparable to Fig. 7 upper right for 
as = 0.8. 



lookback time 
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1:5 
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20% 
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0.0 
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0.16 
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0.12 
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0.33 


0.15 
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6.2 


0.16 


0.22 


0.32 


0.48 


0.26 


346 



Tabic 3 

The number of mergers or large AM events within the previous 6t = 300/i-iMpc 
as a function of lookback time for our as = 1.0 simulation, comparable to Fig. 7 
lower left for (Ts = 0.8. 



lookback time 


1:3 


1:5 


1:10 


20% 


25% 


-^clus 


0.4 


0.06 


0.08 


0.12 


0.13 


0.09 


872 


2.3 


0.08 


0.11 


0.19 


0.21 


0.12 


718 


4.3 


0.10 


0.13 


0.23 


0.33 


0.15 


516 


6.3 


0.14 


0.18 


0.33 


0.57 


0.26 


291 



Table 4 

The number of mergers or large AM events within the previous a^^^ Gyr as a 
function of lookback time for our as = 1.0 simulation, comparable to Fig. 7 lower 
right for (Tg = 0-8. 
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Fig. 8. Top: recent mergers or large AM for a relaxation time of 5t = 600 /i ^Mpc, 
for two values of as- Bottom: the recently merged/large AM fraction for clus- 
ters with M > 3.0 X IO^^^^^Mq, for relaxation times of St = 300^~^Mpc and 
2.5 Gyr. Smooth lines are 1:10, 1:5 and 1:3 predecessor mass ratios, dotted lines are 
Mf/Mi > 1.2, 1.33 (top to bottom). Errors as before. 

and z ~ 0.67 is not uncommon, and sometimes even larger increases occur. 
The errors in these ratios for the relaxation time of 5t — 100/i~^Mpc are 
quite large. For erg = 0.8(1.0), from 44(29)% for the 1:3 mergers to 32(28)% 
for Mf/Mi > 1-2. For the other relaxation times the errors on the ratio are 
the largest for the 1:3 mergers (10-25%) and smallest for the mass jumps with 
Mf/Mi > 1-2 (3-12%). The errors on the other cases fall in between. More 
details can be obtained from the earher plots and the Tables. 



Another question of interest might be: in how many major mergers do three 
or four predecessors have comparable mass? We took one particular definition 
as described in the methods section. These three-body mergers were fairly 
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Fig. 9. Summary of change in recent merger /recent large AM fraction of clusters 
with time: (fraction 6.2 Gyr ago) /(fraction at present) for most, (fraction 6.3 Gyr 
ago) / (fraction 0.4 Gyr ago) for the a^/^ relaxation time case. The various relaxation 
times are listed at the bottom: lOOt, 300t, 600t {5t = 100, 300, 600 /I'^Mpc respec- 
tively), tcross = CL^^"^ Gyr, and 1.0 and 2.5 Gyr. Filled triangles are 1:3 mergers, open 
squares are 1:5 mergers, filled squares are 1:10 mergers, open circles are 20% mass 
jumps, filled circles are 33% mass jumps. The errors are fairly large at the early 
time but one can see that a doubling of the recently merged or large AM frac- 
tion between the present and z ~ 0.8 is not uncommon and sometimes even larger 
increases occur. 

common, the exact numbers depend upon relaxation times considered. For 
the CTg = 0.8 case, about 15% of the 1:5 mergers were three-body for 6t = 
30 h~^Mpc and slightly less than 30% for 6t = 600 h~^Mpc on average, with 
little evidence of growth with time. For erg = 1.0 the fraction of three-body 
mergers to 1:5 major mergers was closer to 10% on average and seemed to grow 
with time for St = 300/i~^Mpc. For 6t = 600 /;,~^Mpc this growth was more 
evident, starting around 12% and growing to around 30%. For fixed relaxation 
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times, 2.5 (1.0) Gyr, the fraction of three-body major mergers to 1:5 major 
mergers showed a definite increase with lookback time. For (jg = 0.8 it started 
at 30% (10%) and grew to 80% (35%) at lookback times of 7 Gyr. For ag = 1.0 
the fraction started at 12% (5%) and grew to 55% (25%) at lookback times 
of 7 Gyr. This means that even though a 1:5 merger might not be considered 
"major" by some definitions, it is very possible that there are two of these 
going on at the same time, making these mergers more energetic. 

There were many fewer four-body mergers. The largest effect was for the 
2.5 Gyr relaxation time, where four-body mergers were about one third of the 

three-body mergers. In the largest case, for cig = 0.8 and 2.5 Gyr relaxation 
time, the total fraction of clusters which had had a recent four-body merger 
only reached 5% or above for lookback times greater than 4 Gyr. 



5 Conclusions 

Clusters of galaxies represent the current endpoint of structure formation. As 
the largest systems which have had time to virialize in a universe with hier- 
archical structure formation, they make excellent laboratories for cosmology, 
large-scale structure and galaxy formation. The formation of galaxy clusters 
via a combination of mergers and accretion of smaller objects is crucial to 
understanding many of their present day properties. 

In this paper we investigated the formation of galaxy clusters in some detail, 
and the major events which define this process, extending earlier studies (2; 3; 
4; 5; 6; 7; 8; 9; 12) discussed in §2. While galaxy cluster formation properties 
can be reliably calculated with the present generation of large, high-resolution 
N-body simulations, it is difficult to find specific numbers in the literature. 
To remedy this we have calculated, for a sample of hundreds of clusters, the 
degree of virialization, formation time, and for two values of the clustering 
amplitude ag, the fraction of disturbed galaxy clusters for many definitions of 
disturbed. 

We began by showing the time histories of a few clusters to illustrate the 
typical formation pattern: periods of smooth accretion punctuated by large 
increases in mass. We then turned to various characterizations of this process, 
applied to our statistical sample as a whole. 

We first calculated the ratio of kinetic to potential energy. At early times 
clusters are hotter than the vacuum virial relation (2KE=PE) would predict 
due to continuous infall of material. With the onset of A domination and the 
cessation of structure growth the excess kinetic energy drops and their mass 
accretion rate slows. 
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We then turned to formation time definitions in the hterature; while no sin- 
gle number can capture the complexity of a cluster merging tree, the concept 
of a formation time encodes much useful information. We compared difi:erent 
formation time definitions which are relevant for different types of observa- 
tions. For instance a recent large mass jump may not be of interest if one only 
wants to know when the cluster first was detectable by SZ decrement (which 
presumably depends more on the depth of the potential well), but may be rel- 
evant for studies which rely on a cluster being old enough to be dynamically 
relaxed. The different formation times we considered are correlated, with large 
scatter. Formation times relying upon a smooth parameterized fit to galaxy 
cluster histories were most correlated with each other, but the results were 
sensitive to the fit methodology, due to the large amount of recent sporadic 
mass gain. We compared the distributions of the different formation times for 
our whole cluster sample and for our more massive clusters, again finding some 
correlation but also significant scatter. 



We then turned from smooth parameterizations to the characteristic abrupt 
jumps in mass over time for galaxy clusters. On average, galaxy clusters have 
had at least 4 large mass jumps since ^ ~ 2, this number increases with cluster 
mass; about half of the clusters get at least half of their mass in these jump 
events. 



The above measurements were for our (Xg = 0.8 sample of 574 clusters. The 
second part of the paper reported the time dependence of the fraction of 
clusters which have had a recent mass jump or major merger, using in addition 
a (jg = 1.0 sample of 909 clusters. These fractions should be of use for estimates 
of the number of "relaxed" clusters available in surveys (given a relaxation time 
for the phenomena of interest) and for helping to constrain relaxation times for 
phenomena associated with mergers or mass jumps where the phenomenon's 
occurrence fraction has been measured but not its relaxation time. In Figs. 7-9 
and tables 1-4 we give the merger fractions for a number of different situations. 
These fractions also serve to quantify the previously seen trend of more recent 
mergers (fractionally) at higher redshift. Comparing z = with z = 0.67 
(about 6 Gyr ago), the fraction of clusters which had a recent merger or mass 
jump is doubled by almost any definition and for some definitions the increase 
is even larger. 



The simulations used here were performed on the IBM-SP at NERSC. JDC 
thanks T. Abel, S. Allen, G. Bryan, R. Gal, J. Hennawi, R. Kneissl, A. 
Kravtsov and P. Ricker for helpful discussions and was supported in part 
by NSF-AST-0205935. MJW was supported in part by NASA and the NSF. 
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